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We propose an alternative formulation of structure functions for the velocity field in fully de- 
veloped turbulence. Instead of averaging moments of the velocity differences as a function of the 
distance, we suggest to average moments of the distances as a function of the velocity difference. 
This is like an "inverted" structure function, with a different statistics. On the basis of shell model 
calculations we obtain a new multiscaling spectrum. 



ON 
ON 
ON 

< 

m 
> 

(N : 
o ■ 

o : 

ON . 
ON ■ 

"c : 

>v 
■ 

i 1 

o : 
^ . 

o ■ 



- 1—1 

X 



The understanding of intermittency effects in fully de- 
veloped turbulence and the associated multiscaling spec- 
trum of exponents, is probably the most fundamental 
open problem in turbulence research Q . The traditional 
way of describing this is, as already suggested by Kol- 
mogorov [Q, to consider the velocity difference between 
two points of the turbulent state, raise this difference 
to the moment q, and then study the variation with re- 
spect to the distance between the two points. To im- 
prove the statistics, the moments are averaged in space 
(and maybe time) and one obtains the well known struc- 
ture functions where the corresponding scaling exponents 
are called structure function exponents Q. During the 
last decades it has become clear both from many experi- 
mental , numerical || and theoretical considerations 
j^D , that this set of exponents is very non-trivial, defin- 
ing an infinity of independent exponents leading to at 
"curved" variation of the scaling exponent with the mo- 
ment. Notable is also the recent fundamental advances 
in obtaining the multiscaling spectrum analytically for 
a passive scalar advection in a spatially correlated, but 
temporally non-correlated, velocity fields, the socallcd 
Kraichnan model ||[l0| . 

We propose simply to "invert" the structure function 
equation, and consider instead averaged moments of the 
distance between two points, given a velocity difference 
between those points. This leads to an alternative way 
of describing and analyzing a turbulent velocity field (in 
particular when measured experimentally) and one ob- 
tains a new set of exponents that we have not yet been 
able to relate to the traditionally estimated exponents, 
though we suspect that there might be a relation. This 
inversion is inspired by studies in passive scalar advection 
where one often, say for pair particles, considers averages 
of the advection time versus the distance, instead of av- 
erages of the distance versus time Jll]-[l4|] . To a start let 
us introduce the well known structure functions for the 
velocity field u(x, t) of a fully developed turbulent state, 
obtained either from the Navier-Stokes equations or from 
measurements 



Alternatively, we now consider the following quantities, 
which is denoted the distance structure functions 



< Au x (£) q >~ f« 
where the difference is defined as 
Au x (£) = u(x + r) - u(x) 



(1) 



(2) 



< £(Au x ) q >~ |Ai 



(3) 



where the difference Au x is again defined as in Eq. (g) 
and £{Au x ) is understood as the minimal distance in r, 
measured from x, for which the velocity difference ex- 
ceeds the value Au x . In other words, we fix a certain 
set of values of the velocity difference Ait x . Starting out 
from the point x, we monitor the distances £(Au x ) where 
the velocity differences are equal to the prescribed values. 
Performing an average over space (and maybe time) the 
distance structure functions Eq. (||) are obtained. By 
assuming self-similarity of the small scale velocity differ- 
ences, one expects a trivial set of exponents 5 q where the 
variation with the moment q is determined by one expo- 
nent. Say, in the standard Kolmogorov theory we know 
that the velocity differences behave as Au ~ f 1 / 3 , forget- 
ting for a moment the averaging brackets. Inverting this 
equations, we of course obtain £ ~ Au 3 and would ex- 
pect a trivial relation 6 q = 3q. In case of an intermittent 
and singular velocity field without self-similarity of the 
small scale velocity differences (see |l5)), this would be 
completely different and the averaging brackets will be 
crucial, relating to the statistics of the varying quantity 
that is averaged. We will show, based on shell model cal- 
culations, that in turbulence there exists a new spectrum 
aears not to be trivially related to the spec- 
,17|| . Let us for a moment reflect on the case 



S q , that a 
trum C, q ]l 



The average in Eq.(Q) is over space (and maybe time) . 
We have assumed full isotropy of the velocity field. The 
set of exponents C, q forms a multiscaling spectrum M . 



q = 1. Using the standard value £i ~ 0.38 — 0.40, the 
simple inversion gives 5\ ~ 2.5. Our calculations indicate 
that this value is not obtained in a turbulent model field. 
Instead we find a value Si ~ 2.0 — 2.1. Another way to do 
the comparison is to aim at velocity exponent 1 and find 
the corresponding moment q, i.e. < £ g >~ Au 1 . We ob- 
tain q ~ 0.45, again different from 0.40. These differences 
are of course attributed to the very different statistics, 
i.e. whether the velocity differences or the corresponding 
distances are averaged. Also, we obtain strong intermit- 
tency corrections in the sense that the value of the 8'th 
moment is Cs ~ 12.9, i.e. much smaller than 8£i ~ 16.3. 

In order to apply this scheme in a direct calculation 
we employ the Gledzer-Ohkitani-Yamada, GOY, shell 
model |l^,19| which has be intensively studied over the 
last years p0|-|26]|. This model is a rough approximation 
to the Navier-Stokes equations and is formulated on a 
discrete set of k- values, fc„ = r n . We use the standard 
value r = 2. In term of a complex Fourier mode, u n , of 
the velocity field the model reads 



1 



( i 'hi * * * . 

{— +VK n ) U n — I K n (a n U n+1 U n+2 + "7f w n-l u ri+l ' 



(4) 



with boundary conditions 61 = = c\ — C2 = oln—i = 
a at = 0. / is an external, constant forcing, here on the 
forth mode. 
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FIG. 1. a): The velocity structure function of order one. 
The line has a slope of 0.39. b): The distance structure func- 
tion of order one. The line has a slope of 2.02. Note the inner 
cut-off related to the dissipative cut-off in a), and the outer 
cut-off given by velocity of the forcing scale, c): The distance 
structure function of order 8. The exponent is Ss ~ 13.1. The 
"raggedness" is due to discretization of the varying length 
scale. 

The coefficients of the non-linear terms must follow the 
relation a n + b n+ i + c n+ 2 = in order to satisfy the con- 
servation of energy, E — J2n |u„| 2 , when / = v = 0. 
The constraints still leave a free parameter e so that 
one can set a n = 1, b n+ i = -e, c n+2 = -(1 - e) 
p3 . As observed by Kadanoff , one obtains the canonical 



value e = 1/2 , if helicity conservation is also demanded 
p3f . The set (||) of N coupled ordinary differential equa- 
tions can be numerically integrated by standard tech- 
niques. We have used standard parameters in this paper 
N = 27, v = 1(T 9 , k = 0.05, / = 5 • 1CT 3 . 
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FIG. 2. a): The exponents £ g for the velocity structure 
functions, with selected error bars. The line corresponds to 
Kolmogorov theory, b): The exponents S q for the distance 
structure functions. The line is adjusted to pass through the 
value of the first order exponent (l,<5i). 

The GOY model is defined in fc-space but our formal- 
ism is written in direct space and we therefore apply a 
sort of inverse Fourier transfor m p7| . Here we employ 
an idea proposed by Vulpiani p8| , |29| ] and write the three- 
dimensional velocity field in the following way 



N 



U(X, t) = Cn [Un {t)^ n ' X + C. c] 



The wavevectors are defined by 



(5) 



(6) 



where e„ is a unit vector in a random direction, for each 
shell n. Also c„ are unit vectors in random directions. 
One can easily ensure that the velocity field is incom- 
pressible, div u = 0, by the following constraint p8| 



e„ =0 Vn 



(7) 



Note, that this condition could be relaxed to XmLi c « ' 
e„ = 0. In our numerical computations we consider the 
vectors c n and e„ quenched in time but nevertheless av- 
erage over many different realizations of these; i.e. one, 



2 



or several, specific measurements of the distance struc- 
ture functions are performed with one realization of the 
vectors. After that a new realization of e„,c„ is applied 
in order to perform a good statistical average. 
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FIG. 3. Probability distribution 

functions P(£(Au*)\Au x ) for (a) the velocity Au x = 0.0027, 
which is close to the dissipative length scale (see Fig.l), and 
for (b) the velocity Au x = 0.26, close to the velocity of the 
outer cut-off. 

Equipped with a real space time dependent velocity 
field we start out with a test of this field by computing 
the standard velocity structure functions, given by Eq. 
(0) ©■ Indeed, the field exhibits nice scaling invariance 
as shown in Fig. la, where the first order velocity struc- 
ture function is presented. We have extracted all the 
exponents with moment up to g=10 and the correspond- 
ing results are shown in Fig. 2a (and for completeness also 
in Table 1). These results agree with the exponents ob- 
tained b y n umerical computations of the GOY-model in 
/c-space pql , i.e. without performing the transformation 
to real space. In the averaging, we have assumed isotropy 
and for practical convenience, the distance is varied only 
along the three coordinate axes. Having checked this we 
proceed to extract the distance structure functions, Eq. 
(||). Practically, both the distance and velocity differ- 
ences are discretized as £ = A^ and Au = Aj. In the 
present calculations the value Xd = A u = 1.02 is cho- 
sen. As the starting point we set x = and vary again 
along the coordinate axes. For a fixed value of Aito, t is 
increased until for the first time the velocity difference ex- 
ceeds this fixed value: this defines £(Aqu). Then Auo is 
increased by one more step and so on. Fig. lb presents the 
scaling of the first order distance structure function and 
the corresponding exponent 8i is estimated to a rather 



good precision, 6i — 2.02 ± 0.05, with a scaling regime 
of 2 decades on the A u axes and 4-5 decades on the £ 
axes. Note, the cut-off at low values of Aou. This cut-off 
is related, both for values of velocity and distance, to the 
dissipative cut-off of the standard structure function, see 
Fig. la. The cut-off at large values of Aqu is related to 
the velocity at the forcing scale. In all the presented cal- 
culations we have averaged over 24630 situations. Fig. lc 
presents the distance structure function of order q = 8, 
resulting in an exponent Sg = 12.9 ±0.5. The graph is 
"rough" due to the binning of Aou and due to the high 
moment. Fig. 2b shows the multiscaling spectrum of 6 q . 
We have included a straight line through the point (1, Si) 
in order to show the curved nature of the spectrum. For 
completeness, Table 1 also displays the measured scaling 
exponents 5 q . 

It is well known, that one can improved the scaling sig- 
nificantly using the technique of extended self similarity 
(ESS) |3l| where one moment of a given variable is varied 
against another moment. In the present case this means 
a graph of one distance structure function < £{A ) q > 
versus another < £(Ao) q > for two different moments 
q, q' , and this results in ESS plots which spans over three 
times as long a regime as compared to traditional ESS 
plots where the quantities are < Au x (£) q > are applied 
(the large regime is of course due to the Kolmogorov | 
exponent relation). Applying ESS we have obtained the 
exponents 8 q in an independent way and the results agree 
well with the values listed in Table 1. This property of 
a much larger scaling regime of the ESS plots could be 
one of the advantages of the presented formalism. Details 
will be given in a forthcoming publication. 

To obtain a better understanding of the obtained re- 
sults we need to consider the statistics of £(Au x ) q in the 
following way 

< £(Au x ) q >~ / £(Au x ) q P{£(Au x )\Au x ) d£ (8) 



where we have introduced the conditional probability dis- 
tribution function P(£(Au x )\Au x ). This measures the 
probability of a distance £ given the velocity difference. 
We show this PDF for two different values of the veloc- 
ity difference in Fig. 3 on linear scales. In both cases, 
the distributions are clearly non-Gaussian with long ex- 
ponential (or in fact stretched exponential) tails, as ex- 
pected in intermittent systems. The surprising difference 
to the standard PDF's for velocity differences is, that it 
does not tend towards a Gaussian for large scales. We 
would have expected that. We have not been able to 
relate this PDF, P(£(Au x )\Au x ), to the "usual" PDF, 
P(Au x \£); these two PDF's measure simply very differ- 
ent things. 

In conclusion, we have introduced the distance struc- 
ture functions defined for a velocity field in fully devel- 
oped turbulence. The corresponding multiscaling spec- 
trum appears not to be related to the well known spec- 
trum for velocity structure functions. The distance struc- 
ture function could be very relevant for experimental ve- 
locity data measured in one point fl7| | . Here one typically 
applies the Taylor hypothesis in order to relate a tempo- 



3 



ral segment to a spatial segment. For this type of time 
series, the distance structure functions should be easily 
extracted. 
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TABLE I. Values of the scaling exponents for velocity 
structure functions £ 9 and distance structure functions 8 q with 
selected error bars. 



q 


1 


2 


3 


4 


5 


6 


7 


8 


5 q 


0.39 
2.04 


0.73(2) 
3.70(5) 


1.0 
5.4 


1.28(5) 
7.0(2) 


1.53 
8.53 


1.77(6) 
10.0(4) 


2.0 
11.7 


2.20(8) 
12.9(6) 
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